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ABSTRACT 

We present Hubble Space Telescope (HST) and ground-based observations of 
a pair of galaxies at redshift 2.38, which are collectively known as 2142—4420 Bl 
(Francis et al. 1996). The two galaxies are both luminous extremely red objects 
(EROs), separated by 0.8". They are embedded within a 100 kpc scale diffuse 
Lya nebula (or blob) of luminosity ~ 10*^^ erg s~^. 

The radial profiles and colors of both red objects are most naturally explained 
if they are young elliptical galaxies: the most distant yet found. It is not, however, 
possible to rule out a model in which they are abnormally compact, extremely 
dusty starbursting disk galaxies. If they are elliptical galaxies, their stellar pop- 
ulations have inferred masses of ~ IO^^Mq and ages of ~ 7 x 10* years. Both 
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galaxies have color gradients: their centers are significantly bluer than their outer 
regions. The surface brightness of both galaxies is roughly an order of magnitude 
greater than would be predicted by the Kormendy relation. A chain of diffuse 
star formation extending 1" from the galaxies may be evidence that they are 
interacting or merging. 

The Lya nebula surrounding the galaxies shows apparent velocity substruc- 
ture of amplitude ~ 700 km s~^. We propose that the Lya emission from this 
nebula may be produced by fast shocks, powered either by a galactic superwind 
or by the release of gravitational potential energy. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: interac- 
tion — galaxies: individual (2142-4420 Bl) 

1. Introduction 

Two of the most enigmatic types of high redshift galaxy are the extremely red objects 
(EROs, eg. Hu & Ridgway 1994; Thompson et al. 1999), and the Lya blobs (eg. Lowenthal 
et al. 1991; Francis et al. 1996; Keel et al. 1999; Steidel et al. 2000; Roche, Lowenthal & 
Woodgate 2000; Kobulnicky & Koo 2000). 

EROs are defined as having extremely red observed- frame optical/near-IR colors, and 
may be dusty starburst galaxies or AGN (eg. Small et al. 1999; Hughes et al. 1998), or 
spheroidal galaxies (eg. Cimatti et al. 1999; Dunlop et al. 1996; Soifer et al. 1999; Moriondo, 
Cimatti & Daddi 2000). 

The Lya blobs are characterized by 100 kpc scale low surface brightness Lya nebulae, 
with total Lya luminosities > IC^'^crg s~^. All Lya blobs identified to date appear to lie 
in proto-cluster environments. It has been proposed that these blobs are a form of AGN, 
that they are associated with cooling flows (Steidel et al. 2000; Haiman, Spaans & Quataert 
2000) or that they are generated by galaxy superwinds (Taniguchi & Shioya 2000). 

Can a galaxy be both an ERG and a Lya blob? Many high redshift radio galaxies 
contain ERG nuclei embedded within dense gaseous regions, which are often associated with 
extended diffuse Lya emission (eg. McCarthy 1993; Pentericci et al. 1997; Ivison et al. 2000; 
Garilh et al. 1997; Binette et al. 2000; BickneU et al. 2000). 

In this paper we present a detailed study of a radio quiet Lya blob that contains two 
compact ERGs. Francis et al. (1996) identifled 2142—4420 Bl, which is a luminous Lya 
and C IV emitting source at redshift z = 2.38. The enormous Lya luminosity of this source 
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comes from an extended ~ 100 kpc low surface brightness nebula. Ground-based near-IR 
observations showed an unresolved, luminous, extremely red source lying within this blob. 
The whole system appears to lie within a cluster of Lyman-limit QSO absorption-line systems 
(Francis, Wilson & Woodgatc 2000). The optical properties of this source closely resemble 
those of many high redshift radio galaxies, but Francis et al. (1996) obtained an upper limit 
on its radio flux of 0.27 mJy at 1344 and 2378 MHz. 

We describe our observations in Section 2, and show the results in Section 3. We then 
discuss the nature of the ERO (Section 4) and of the Lja blob surrounding it (Section 5). In 
Section 6 wc conclude that this object is probably a pair of young merging elliptical galaxies, 
though a compact dusty starburst model is hard to completely exclude. We further conclude 
that the Lya nebula surrounding them is either photoionized by an AGN or shock excited 
by the energy of gravitational infall or of a superwind. 

We assume a fiat universe with Hq = 70km s~^Mpc~^, = 0.3 and VIa = 0.7. At 
z — 2.38, given this cosmology and redshift, one arcsecond corresponds to 8.0 proper Kpc, 
and the luminosity distance is 18.8 Gpc. 

2. Observations 

Our previous observations of this object can be found in Francis & Hcwctt (1993); Fran- 
cis et al. (1996); Francis, Woodgate & Banks (1997) and Francis, Wilson & Woodgate (2000). 
The observation log is in Table 1. HST imaging was used to determine the morphology and 
colors of Bl. Ground-based imaging measured its integrated colors and diffuse Lya emission. 
Shtless spectroscopy and Fabry-Perot imaging was used to constrain the velocity structure. 

In addition to our own observations, we searched the ROSAT All Sky Survey (Vogcs et 
al. 1999) for data on this object. There is only a 308 sec exposure toward the field, giving 
a 90% upper hmit X-ray count rate of 0.0075 counts sec"^. Assuming a Galactic column 
density of 2.65 x 10^*^ cm"^ along this fine of sight (which is the weighted mean of the four 
nearest 21 cm measurements) and a power law spectrum with an energy index of 1.0, we 
obtain a flux /e < 2.6 x 10~^ photons cm~^ sec"^ keV~^ where E is the energy in keV. 

2.1. WFPC2 Imaging 

The cluster region was observed for 38 orbits with the Wide Field and Planetary Camera 
II (WFPC2, Trauger et al. 1994) on the Hubble Space Telescope. Galaxy Bl was placed 
on the WF3 chip. Exposures were dithered on a sub-pixel grid, to allow reconstruction of a 
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Tablc 1. Observations 



Instrument 


Filter 


Rest- frame z — 2.38 
wavelength (nm) 


Exposure 
(s) 




Date 




Notes 


HST/WFPC2 


F410M 


121 ±2.1 


67200 


1999 


Aug 


27- 


-Sep 8 




HST/WFPC2 


F450W 


131 ± 13 


19400 


1999 


Aug 


31- 


-Sep 2 




HST/WFPC2 


F702W 


207 ±21 


14400 


1999 


Aug 


23 






HST/NICMOS 


FllOW 


325 ± 87 


5120 


1998 


Jun 


27 






HST/NICMOS 


F160W 


471 ± 54 


10112 


1998 


Oct 


30 




a 


HST/NICMOS 


F164N 


487 ±2.5 


10236 


1997 


Oct 


30 




a 


HST/STIS 


Clear 




3120 


1998 


Oct 


25- 


-27 


b 


HST/STIS 


G430L 




17000 


1998 


Oct 


25- 


-27 


c 


NTT/SOFI 


J 


369 ± 43 


3540 


1999 


Aug 


16 






NTT/SOFI 


H 


489 ± 44 


4740 


1999 


Aug 


16 






NTT/SOFI 


Ks 


640 ± 40 


3480 


1999 


Aug 


16 






CTIO/Mosaic 


4107A 


121.5 ±0.8 


18000 


1999 


Aug 


7 






CTIO/Mosaic 


U 


108 ± 10 


9000 


1999 


Aug 


8 






CTIO/Mosaic 


B 


130 ± 14 


7200 


1999 


Aug 


7- 


-8 




CTIO/Mosaic 


V 


163 ± 13 


7800 


1999 


Aug 


7- 


-8 




CTIO/Mosaic 


R 


207 ± 32 


2400 


1999 


Aug 


8 






CTIO/Mosaic 


I 


260 ± 35 


4800 


1999 


Aug 


8 






AAT/TTF 


4110A 


121.6 ±0.6 


5400 


1999 


Sep 


11 




d 


AAT/TTF 


4114A? 


121.7 ±0.6 


5400 


1999 


Sep 


11 




d 


AAT/TTF 


4118A 


4121.8 ±0.6 


9000 


1999 


Sep 


12 




d 



*Not taken during a NIC3 focus campaign. Image quality degraded, and image partially 
vignetted 

'^Taken to allow registration of the slitless spectroscopy 
'^Slitless spectroscopy 
^Wavelength calibration uncertain 
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better sampled image. Broad-band images were obtained through the F450W and F702W 
filters. Relatively narrow-band images were taken through the F410M filter, which matches 
the wavelength of Lya at the cluster wavelength. The F410M images suffer from a very low 
surface brightness, time-varying mottling. We hypothesise that this is caused by scattered 
Earth-light. 

The images were pipeline processed, then cosmic rays were removed and the frames 
co-added using the drizzle algorithm (Fruchter & Hook 1998) to minimise undersampling. 

2.2. STIS Slitless Spectroscopy 

The Space Telescope Imaging Spectrograph (STIS, Woodgate et al. 1998) observations 
(pairs of slitless spectra and images) were taken at two different orientation angles separated 
by 11.5 deg, in an attempt to separate slitless spectra from objects aligned along the same 
row for a given angle. For both the direct clear-aperture and the dispersed images, the CCD 
was binned 2x2, producing a plate scale of 0.1 arcsec/pixel, and spectral resolution ~ 11 
A for point sources. There were 20 exposures each for the direct and dispersed images, with 
dithering of ~ 1.4 arcsec between each exposure in order to minimize the effects of cosmic 
rays and hot pixels. The direct and dispersed images were combined using the STIS GTO 
team software CALSTIS (Lindler 1998; Gardner et al. 1998). The spectra were extracted 
using both the STIS_extract (Chen, Lanzetta & Pascarelle 1999) and SLWIDGET (Lindler 
2000) software. Results from both were consistent. 

2.3. NICMOS Imaging 

The Near Infrared Camera and Multi Object Spectrometer (NICMOS, Thompson et 
al. 1998) observations were made with Camera 3, through the FllOW, F160W and F164N 
filters. All exposures were dithered for each filter over 4 (or 8 in the case of F160W) positions 
spanning a 2 x 2 arcsec box, to minimize the effects of bad pixels. The F160W and F164N 
observations are affected by vignetting from the field offset mirror, and were taken before 
the focus became optimal. The NICMOS data were reduced using the same methods as for 
SOFl (Section 2.4). The data were also reduced using NICRED (McLeod 1997); there is little 
difference between final images reduced by the two methods. A comparison of broadband 
SOFI and NICMOS data show the photometry to be consistent between the two, with the 
exception of a vignetted band along one end of the chip in the F160W and F164W images. 
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2.4. NTT Imaging 

Near-IR photometry was obtained with the SOFI camera on the ESO New Technology 
Telescope (NTT, Lidman, Cuby & Vanzi 2000) . Conditions were photometric, but seeing 
was ~ 1.1" and the telescope was subject to considerable wind shake. The data were reduced 
using a modified version of IRAF^^ scripts written by Peter McGregor. NICMOS standards 
(Persson et al. 1998) were used. 

2.5. Cerro Tololo 4m Imaging 

We obtained broad- and narrow-band (Lya) imaging of this field with the MOSAIC 
camera (MuUer et al. 1998), at the prime focus of the Cerro Tololo Inter- American Observa- 
tory (CTIO) 4m Blanco Telescope. The narrow-band images were taken through a specially 
purchased filter, which working in the f/2.8 beam has a bandpass of 54A, centered at 4107A. 
Conditions were mostly photometric, with seeing ~ 1.2". 

2.6. Anglo-Australian Telescope Narrow-band Imaging and Spectroscopy 

Our Anglo- Australian Telescope imaging was obtained with the Taurus Tunable Filter 
(TTF, Bland-Hawthorn & Heath- Jones 1998), which is a Fabry- Perot etalon system. Con- 
ditions were photometric on the first night, and the typical image quality was ~ 1.6". A 
MIT-Lincoln Labs CCD was used behind the TTF. 

The TTF was used with a full-width at half maximum (FWHM) spectral resolution of 
between 3.35 and 4.4A. Images were obtained at central wavelengths of 4110A, 4114A, and 
4118A. 

An electronics fault was subsequently found to induce significant wavelength drifts on 
timescales of a few hours. The 4110 and 4118A observations were taken within four hours 
of a cafibration, but the 4114A observations were taken between four and eight hours after 
the calibration. Thus the wavelength calibration of all three images are somewhat uncertain, 
particularly that of the 4114A image. 

We obtained a long-slit spectrum of Bl with the Low Dispersion Survey Spectrograph 



^^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation. 
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(LDSS, Wynne & Worswick 1988). These observations are described in Francis, Woodgate 
& Danks (1997) 

3. Results 

Our images of Bl are shown in Figures 1 & 2. The appearance of Bl varies dramatically 
as a function of wavelength and resolution. In the observed-frame red and near-IR, it breaks 
up into two compact components: Bla and Bib (Fig 2), separated by 0.8". In this respect 
it is similar to another strong, high redshift Lycc source, the Coup-Fouree galaxy (Roche, 
Lowenthal & Woodgate 2000). A third continuum source, Blc, appears as chain of diffuse 
knots extending about 1" to the south from Bla and Bib. 

In Lya, the picture is different again. When viewed with WFPC2, we see a knot of 
Lya emission from Bib, but no emission from Bla. Diffuse patches of low surface brightness 
Lya emission are found ~ 1" to the south and north-east of Bib. The ground-based CTIO 
image has much greater sensitivity to low surface brightness Lya emission, and confirms 
that it extends at least 10" north-east of Bla and Bib, as reported by Francis et al. (1996) 
and Francis, Woodgate & Danks (1997). The non-detection of this emission in the WFPC2 
image indicates that it is truly diffuse, and not coming from a series of compact sources 
separated by less than the ground-based resolution limit (as was the case for the Lya nebula 
observed by Pascarelle et al. (1996)). 

The LycK contribution to the blue continuum images, and the blue continuum contri- 
bution to the narrow-band Lya images, were removed. Bl has a total Lya flux of about 
1.5 X 10~^^erg cm~^s~^: a more exact number is hard to calculate as the diffuse flux fades 
gradually into the sky (this number is for a rectangular aperture 10" long and 5" wide, aligned 
along the major axis of the Lya flux). Bib has a Lya flux of ~ 2.0 x 10~^^erg cm~^s~^, 
which is only ~ 13% of the total. The remaining fiux is diffuse: the average Lja surface 
brightness of the region within 2" of Bib is ~ 4 x lO^^^erg cm~^s~^arcsec~^. 4" further east, 
the average surface brightness drops to ~ 10~^^erg cm~^s~^arcsec~^. 

3.1. Fluxes 

Continuum fluxes (Table 3) were calculated for the three components of Bl. Small 
circular apertures were used: 0.3" in radius for Bla and Bib, and 0.65" for Blc. These 
small apertures will miss much of the flux of these components, but were chosen to avoid 
contamination from the other components. No aperture corrections were attempted. A total 
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Fig. 1. — Images of a 20" by 15" region around the centre of Bl. The bottom left image is 
the CTIO Lya image. The F410M and F450W images have been smoothed with a Gaussian 
of 0" = 2 pixels (0.1"). The CTIO image has been rotated to match the orientation of the 
WFPC2 images. The continuum has been removed from the Lya images, and the Lya 
component has been removed from the F450W image. 
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Fig. 2. — Close-up of the continuum sources in Bl. The F410M and F450W images have 
been smoothed with a Gaussian oi a — 1 pixels (0.1"). The NICMOS images have been 
rotated to match the orientation of the WFPC2 images. The orientation matches that of 
Fig 1. 
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flux for Bl was also measured, using a 2"radius aperture. 

Error estimates were made using apertures placed randomly on sky regions. A copy of 
the image was made with all the sources detected by the SExtractor (Bertin & Arnouts 1996) 
package masked out. We randomly placed 1000 apertures of identical radius on the image. 
Any aperture that fell on a masked region was rejected, and the procedure continued until 
1000 clean measurements had been obtained. The standard deviation of these measurements 
is our sky error estimate. A Poisson noise model was combined with the sky error estimate to 
calculate object flux errors. The Lja contribution to the F450W flux has been removed, but 
the probable emission-line contributions to FllOW, F160W and Kg have not been removed: 
these contributions are discussed in Section 3.1.1. Finally, the fluxes were corrected for the 
local Milky Way extinction {E{B — V) = 0.019), as estimated by Schlegel, Finkbeiner & 
Davis (1998). The extinction curve of CardeUi, Clayton & Mathis (1989) was used. 

Integrated near-IR magnitudes for Bl were measured from the NTT image, using SEx- 
tractor with 4.2" radius circular apertures. We measure J = 21.71 ± 0.12, H — 20.31 ± 0.08 
and K, = 19.47 ±0.10. 

The measured coordinates of the components of Bl are shown in Table 2. 

3.1.1. Emission-line Contamination 

The near-IR fllters are subject to possible emission-line contamination at this redshift, 
z=2.38. [O II] (3727 A) is redshifted into the J-band, H/3 and the [O III] (4959 and 5007 A) 
doublet are shifted into the i7-band, and Ha and [N II] (6583 A) are shifted into the K 
band. These lines could in principle be very strong (eg. Bales et al. 1993). In this section, 
we estimate their strength, and conclude that it is probably small. 

As described in Francis et al. (1996), we imaged the field with a narrow-band filter of 



Table 2. Bl Component Coordinates 



Component 


Coordinates (J2000) 


Bla 


21:42:27.45-44:20:28.69 


Bib 


21:42:27.51-44:20:28.99 


Blc 


21:42:27.46-44:20:30.2 
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wavelength coverage 2.238 ± 0.024//m, which includes both Ha and [N II]. Using a 5" radius 
aperture, wc determine that these lines contribute 21 ± 10% of the integrated flux of Bl. 
The F164N NICMOS image covers the wavelength of redshifted H/5. No significant narrow- 
band excess was detected: we can place a 3cr upper limit on the fraction of the integrated 
i7-band light coming from H/3 in Bl of 3% . 

We have no direct measurements of the other line fluxes, so we estimate their importance 
by making assumptions about their ratios with measured Balmer lines. Typical line ratios 
were taken from Osterbrock & Martel (1993), Kennicutt (1992) and Teplitz et al. (2000), 
for a variety of assumed chemical compositions and ionization sources, and including both 
narrow- and broad-line AGN ratios. 

We conclude that line emission probably contributes only ~ 20% at most of the total 
continuum flux of Bl in H and Kg. In J, the contribution is probably ~ 5% . Furthermore, 
only about 10% of the line flux will likely come from within the photometric aperture of Bib, 
and ~ 20% from within the aperture of the Blc. The remainder of the optical emission hues 
will come from the region of the diffuse blob. As Bla shows no Lya emission, its rest-frame 
optical line flux is also probably small. 

The line corrections will thus probably be minimal. We have not, therefore, applied any 
emission-line corrections to the near-IR fluxes. 



3.2. Velocity Structure 

The Lya line is generally extremely optically thick and often strongly self absorbed, 
making it a very unreliable tracer of dynamics. It is, however, the only sufficiently strong 
line available to us, and can be used to place some limits on the dynamics of the Lya nebula 
surrounding Bl. 

The velocity dispersion of the gas around Bl (measured at the peak of the surface 
brightness, as seen in ground-based images) is ~ 600km s~^ (Francis et al. 1996). Lya 
is redshifted with respect to C IV, indicating probable self absorption. The deep multislit 
spectroscopy of Francis, Woodgate & Banks (1997) marginally spatially resolves this emission 
along the shtlet. Fig 3 shows the spectra extracted from different CCD columns: notice 
the extra component at -l-lOOOkm s~^ in the southernmost spectrum (and another possible 
component at —1000km s~^). 

The diffuse Lya structure differs significantly between the three TTF images. This 
confirms the existence of multiple velocity components separated by ~ 700km s~^. It is 
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Fig. 3. — Longslit Lya spectra of Bl. The slit orientation is shown in Fig 1. The four 
panels are the spectra extracted from the four adjacent CCD columns along which Lya was 
detected: the top panel is the northernmost. Each pixel is 0.83" wide: the seeing was 1.2", so 
the spectra are not independent. Positive velocities represent redshifts. Spectral resolution 
is 700km s~^, so the line core is not resolved in the top three panels. The velocity zero point 
is arbitrary. 
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hard to say more, given the uncertainty in the wavelength cahbration of the TTF data 
(Section 2.6). 

The STIS spectrum detected the Lya flux from Bib. A flux of 1.57±0.40x IQ-^^erg cm-^s-^ 
was measured, consistent with the value derived from the WFPC2 imaging. The hne was not 
signiflcantly spectrally resolved; an upper limit on the velocity dispersion of ~ 1400km s~^ 
can be placed. 

4. The Continuum 
4.1. The Colors 

The spectral energy distributions (SEDs) of the three components of Bl are shown in 
Fig 4. Bla and Bib are quite red, while the SED of Blc is very blue. Blc's colors are similar 
to those of Lyman-break galaxies (eg. Steidel et al. 1996). A power-law fit to its colors gives 
a slope of the form Fx oc A~^'^^, which is close to the median slope of Lyman break galaxies 
(1.5, Adelberger & Steidel 2000). 

Our NTT photometry demonstrates that while Bl is red at wavelengths shortward of 
H, it has a relatively blue H — Kg color. This confirms the result of Francis et al. (1996). 
Bl's SED thus peaks in the if-band. As we discussed in Section 3.1.1, this is unlikely to be 
an artifact of emission-line contamination. 



4. 1.1. Modelling 

We modelled the SEDs of Bl using the 1997 version of the spectral synthesis models of 
Bruzual & Chariot (1993). Models both with continuous uniform star formation and models 
in which all the stars were formed in an instantaneous burst (a simple stellar population) 
were used. All models had Salpeter stellar initial mass functions, with no mass cut-off. Dust 
was modelled using the empirical absorption curves of Calzetti, Kinney & Storchi-Bergmann 
(1994). Our predicted colors are compared with the observations in Figures 5 and 6. 

The red components, Bla and Bib, can be fit by two models: 

• An unreddened burst model of age 750 ± 150 Myr {la limits) and stellar mass ~ 
8 X IO^^Mq. This naturally flts the if-band peak (Fig 6) and the photometry from 
F7021^ thru Kg. It greatly underpredicts the observed F450W flux of Bla. We do not 
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Table 3. Bl Component Observed Flux Densities 



Filter Flux (xlO^^erg cm-V^A ^) 

Bl (total) Bla Bib Blc 



F450W (B) 
F702W (R) 
FllOW (J) 
F160W (H) 



70.2 ± 12.4 
64.2 ±5.1 
47.8 ± 13.6 
120.2 ± 16.4 



11.2 ± 1.3 

12.3 ±0.8 
15.2 ±0.8 
25.7 ±0.9 



< 7.26^ 

10.2 ± 1.0 
11.5 ±0.8 

21.3 ±0.9 



25.7 ±3.4 
15.3 ±2.2 
< 8.35^ 
7.6 ±1.2 



^3(7 upper limit 



10 
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Fig. 4. — The spectral energy distributions of the three components of Bl. Lya emission 
has been subtracted from the bluest point. 
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Fig. 5.— The F450W-F702W and F702W-F160W colors of the three components of Bl. 
Both the continuous star formation models (crosses) and the instantaneous burst models 
(asterisks) have been computed for twelve different ages of the stellar population. From left 
to right, the model ages are 40, 60, 100, 150, 225, 340, 500, 750, 1125, 1700, 2500 and 3800 
Myr. The change in colors caused by dust with a rest-frame V^-band extinction of one 
magnitude is shown. 
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Fig. 6. — The integrated NTT J — H and H — Kg colors of Bl. Model ages and extinction 
arrow as in Fig 5. 
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consider this to be a problem: a star formation rate of only ~ O.lMQyr ^ will produce 
OB stars sufficient to produce the observed blue flux. 

• A dusty starburst model. The best fit is a model with a continuous star formation rate 
of ~ lO^Moyr"^, an extinction of — 2.3 mag, and an age of 500 Myr. A slightly 
worse fit (but still acceptable at the 2a level) is younger (100 Myr), dustier (Ay — 3 
mag) and has a star formation rate of ~ lO^Mgyr"^. 

The dusty starburst model is hard to reconcile with the strong Lya emission of Bib, as 
Lya is resonantly scattered and thus very strongly absorbed by even small amounts of dust. 

Blc is much bluer, and can be fit either by an unreddcncd instantaneous burst model 
of age ~ 100 Myr and mass ~ 3 x IO^Mq, or by a slightly dusty continuous star formation 
model, with a star formation rate of ~ lOMQyr"^ and Ay ~ 1.0. 

4.2. The Radial Profiles 

In this section, we show that Bla and Bib are spatially extended. They are, however, 
very compact objects: too compact to be easily modelled as spiral galaxy disks. They appear 
to show colour gradients: bluer in their central regions. 

In Figures 7 and 8 we compare the radial profiles of Bla and Bib with the relevant 
point spread functions (PSFs) measured from a bright star in the same image. Bla and Bib 
are marginally resolved in the F702W (rest-frame 2100 A) and FllOW (rest-frame 3250 A) 
images. The F160W (rest-frame 4730A) image was taken with NICMOS out of focus, and 
it does not significantly resolve Bla and Bib. The observed PSF does not vary significantly 
across the images, and is consistent with analytic predictions. 

Although Bla and Bib are clearly resolved, both are very compact, with the observed 
surface brightness dropping to half its peak value within about 0.1" (~ 1 kpc). This places 
a strong constraint on the size of the galaxies. Consider a disk galaxy, with a typical radial 
surface brightness profile of the form 

I{r)^IdeM-R/Rd) 

(Binney & Merrifield 1998). The surface brightness of such a disk drops to half its peak 
value at ri/2 = 0.69Rd (Fig 9). Thus our observations require that Rd < 0.15" ~ 1.2 kpc. 
Only ~ 2% of modern galaxy disks are this compact (Kent 1984). Now consider a spheroidal 
galaxy, with a radial surface brightness profile of the form 

7(r) = 7eexp{-7.67[(i?/i?e)'/' - 1]} 
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Fig. 7. — The F702W radial profiles of Bla, Bib and a point source image on the same chip. 
The sohd fine is an analytic fit to the PSF and is the same in all panels. The quadrant of 
each component that faces the other component has been excluded from the plots. 
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Fig. 8. — The FllOW radial profiles of Bla, Bib and a point source from the same image. 
The sohd hne is an analytic fit to the PSF and is the same in all panels. The quadrant of 
each component that faces the other component has been excluded from the plots. 
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(Binney & Merrifield 1998). This profile is much more sharply peaked than the exponential 
disk profile, and its surface brightness drops to half its peak value at ri/2 = 0.13i?e (Fig 9). 
Thus Re < 0.75" ~ 6 kpc. Over 50% of modern elliptical galaxies have effective radii this 
small (Djorgovski & Davis 1987). 

More detailed modelling confirms these results. A wide variety of model galaxy profiles 
were convolved with PSFs and fit by minimisation to the observed radial profiles. The 
modelling confirms that Bla and Bib are significantly extended. Both exponential disk 
models and spheroidal (de Vaucouleurs) models give acceptable fits to the data, as do models 
with a point source embedded in a fainter halo. Pure exponential disk models only give 
acceptable fits if Rd < 0.1" (~ 800 pc), which would place Bla and Bib in the most 
compact 1% of present-day disk galaxies. Spheroidal models give best fits for effective (half- 
light) radii R^, of around 0.5" (~ 4 kpc), which are typical of modern elliptical galaxies. 
Models which combine a nuclear spheroid or point source with a disk can give acceptable 
fits, but only if the disk is unusually small and/or contributes only a small fraction (< 10%) 
of the light. 

The compact radial profiles of Bla and Bib are typical of those of high redshift Lyman- 
break galaxies (Giavahsco, Steidel & Macchetto 1996). The inferred surface brightnesses 
in the rest-frame UV are comparable to those of faint Lyman-break galaxies, but the rest- 
frame optical surface brightnesses are an order of magnitude higher. As Giavalisco, Steidel & 
Macchetto (1996) note, these compact radial profiles are consistent with present-day elliptical 
galaxy profiles or bulges, but are much more compact than present-day disks. 

Curiously, both Bla and Bib appear to be significantly larger at longer wavelengths. 
This can be seen in Figures 7 and 8: despite the wider PSF for the FllOW filter, both 
components are clearly more extended than as seen through the F702W filter. Our modelling 
confirms this result: even after convolution with the relevant PSFs, no single model gives 
a good fit to cither Bla or Bib at both F702W and FllOW. The central ~ 1 kpc of both 
components is ~ 0.3 magnitudes bluer in F702W— FllOW than the outer regions. 

The change in size as a function of wavelength could be explained by a colour gradient 
within a single spheroidal component. It could also be explained by a two component model: 
a compact blue component that dominates at rest-frame 2100A and a more diffuse (but still 
compact) red component that dominates at 3250A and beyond. This red component, if 
disk-hke, must have Rg < 0.2", which would still place it in the most compact 1% of low 
redshift disks. 
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Radius (Arcseconds) 

Fig. 9. — The F702W radial profile of Bla. The dashed curve (left) is a model spheroid with 
Rf. — 0.5". The dotted line (right) is a model disk with — 0.3". Both radii are typical of 
redshift zero galaxies. Neither model has been convolved with the PSF. The solid line is a 
spheroid with Rf. — 0.5", convolved with the PSF and scaled to fit. 
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4.3. Discussion 

What are the continuum components of Bl? Blc, with its blue colors, appears relatively 
straightforward: it is a region of extended moderate star formation, obscured by little dust. 
This star formation may be triggered by the interaction of Bla and Bib. 

The nature of Bla and Bib is less straightforward. We consider four models in turn. 

4.3.1. Elliptical Galaxies 
Could Bla and Bib be young eUiptical galaxies? This hypothesis fits the data well. 

• The radial profiles of Bla and Bib are well fit with de Vaucouleurs profiles, and the 
inferred radii are typical of present-day elliptical galaxies. 

• The red colors are most easily explained by a > IO^^Mq stellar population which 
completed its star formation about 750 Myr before we observe it. It would thus be 
similar to modern E+A galaxies. 

• Bl lies in an overdense region of the early universe (Francis, Wilson & Woodgate 2000), 
which may be the ancestor of a galaxy cluster. 

• The velocity field around Bl, if virial, implies masses of ~ IO^^Mq. 

If Bla and Bib are elliptical galaxies, how do their surface brightnesses and sizes fit on 
the Kormendy relation (Kormendy 1977; Hoessel, Oegerle & Schneider 1987)? The F160W 
filter corresponds quite closely to the rest-frame 5-band. The radial profile fitting suggests 
that both components have effective radii in the observed-frame near-IR of ~ 4kpc. Inte- 
grating a de Vaucouleurs profile over our photometric aperture, we can convert our observed 
F160W magnitudes into rest-frame S-band surface brightnesses at the effective radius. Were 
these components at low redshift, their surface brightnesses at the effective radius would be 
Bqv ~ 19.6 mag arcsec"^. This is roughly an order of magnitude higher than the Kormendy 
relation would predict for an elliptical galaxy of this radius. This is consistent with passive 
evolution from our unreddened burst model (Section 4.1). 

If Bla and Bib are eUiptical galaxies, why do they apparently show color gradients? 
The UV emission from the central ~ 1 kpc could be caused by a star burst, a hole in the dust 
or (for Bib) an AGN. Alternatively, some chemodynamical models predict color gradients 
in young elliptical galaxies (eg. Priaga & Terlevich 1998; Jiminez et al. 1999). 
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4.3.2. Dusty AGN 

A dusty AGN can certainly have colors as red as Bla and Bib (eg. Francis, Whiting 
and Webster 2000). The red components of Bla and Bib are however spatially extended, 
and hence cannot be produced by an AGN. Could the observed color gradient, however, 
be caused by the superposition of a compact blue AGN on a more extended, lower surface 
brightness red galaxy? 

Even in the relatively blue central regions, the continuum slope is redder than oc 
while radio-quiet QSOs have typical continuum slopes of F^, oc (Francis, Whiting and 
Webster 2000). Any central AGN would thus have to be reddened by dust with E{B -V) > 
0.2. This would make it hard for the Lya emission from Bib to escape. In addition, the AGN 
in Bla and Bib would need remarkably equal luminosities and dust extinctions to explain 
the very similar observed colors and radial profiles of the two components. This fine tuning 
seems implausible. Furthermore, the high measured surface brightness was measured in the 
F160W band, to which the central blue component does not significantly contribute. We 
therefore conclude that while AGN may be present, they are unlikely to be responsible for 
the observed colors. Sensitive hard X-ray or far-IR observations will be required to determine 
whether a dusty AGN is present. 



4-3.3. Dusty Disk Galaxies 

The more extended red component of Bla and Bib can be fit by a very compact expo- 
nential disk. Could these two objects then be compact spiral galaxies? The more compact 
bluer central component could be caused by AGNs or nuclear starbursts. 

With sufficient dust, a spiral galaxy can certainly appear as red as these two components 
(Section 4.1.1). Two lines of reasoning however oppose this hypothesis. 

1. Even the extended red components are extremely compact by the standards of modern 
disk galaxies. This could, however, be explained if disk galaxies form from the inside 
outwards. 

2. This model requires that the outer (red) components of Bla and Bib be greatly red- 
dened. This would make it difficult for the observed UV light and Lya emission from 
the bluer central components to escape. 
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4.3.4- Bulges plus Low Surface Brightness Disks 

Could Bla and Bib be the bulges of two disk galaxies? The disks themselves may not 
have formed, or might be too low in surface brightness to be detected by the Hubble Space 
Telescope. 

Two lines of reasoning oppose this hypothesis. Firstly, the red colors and if-band peak 
of Bl are seen in our ground-based photometry with large (~ 5") apertures. Thus any 
extended component must have a similar stellar population to that inferred here. Secondly, 
the inferred masses and sizes for Bla and Bib are far larger than those typical of the bulges 
of disk galaxies. 

5. The LycK Nebula 

Bl has a total Lya luminosity of ~ lO'^^erg s~^, spread over a region at least 30 x 100 
kpc in size. The emission consists of a number of diffuse, discrete components with relative 
velocities of ~ 700km s~^. What can produce such a luminous, diffuse, fast moving Lya 
nebula? 

The bulk motions of the Lya nebula suggest that fast shocks must be present. We will 
show that such shocks are quite capable of producing the observed Lya luminosity. The 
puzzle then becomes one of explaining the origins of the bulk motions. 

It has also been suggested that photoionzation by an AGN, or a cooling flow may power 
such nebulae. We discuss these possibilities, which are hard to exclude. 

5.1. Physical Pcirameters of the Nebula 

If we assume that the Lya velocities are representative of the gas, and not an artifact 
of the high optical depth in Lya, then the dynamical timescale (crossing time) of the nebula 
is ~ 10^ years. This is interestingly similar to the inferred age of the stellar population 
(Section 4.1). 

What is the density of the nebula? The sight-line to background QSO 2139—4434 passes 
20" from Bl (though it does not pass through the Lya nebula). The QSO spectrum shows 
a strong absorption-line system at the redshift of Bl, with a column density Nh ~ 10^^ 
(Francis, Wilson & Woodgate 2000). If we assume that the column density through the 
Lya nebula is at least as great as this, and that the nebula is ~ 20 kpc thick (ie. as thick 
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as it is wide) along the line of sight (both big assumptions), then its density must be at 
least 10~^cm^'^. If, however, the density were this low, the gas would be fully ionized by the 
metagalactic UV background at this redshift. If the nebula is ~ 20 kpc thick, it must have 
a density of > 10~^ for its recombination rate to balance the photoionization from the UV 
background (Francis, Wilson & Woodgate 2000). The diffuse star formation of Blc, taking 
place within part of the Lya nebula, also suggests that densities may be higher. In reality, 
of course, the cloud may be highly inhomogeneous, with dense neutral 'bullets' embedded in 
low density ionized gas. 

If we take this lower limit on the density of the nebula {Nh > 10~^cm~^), and assume 
that the nebula is 100 kpc long, 30 kpc wide and 10 kpc thick, then its hydrogen mass is 
~ 10^°M(^. 



5.2. Shock Models 

Whenever there are blobs of gas moving at supersonic relative speeds, as are observed in 
the nebula, shocks are inevitable. We thus know that shocks are present. Could the emission 
from these shocks produce the observed Lya emission? 

5.2.1. Modelling 

We used the fast shock models computed with the MAPPINGS III code to calculate the 
line emission from the shocks. MAPPINGS III is an updated version of the MAPPINGS II 
code (Dopita & Sutherland 1996) and uses improved radiation continuum and hydrogenic 
calculations to calculate the line and continuum radiation from shocks. The new version 
of the code enabled the calculation of a grid of shocks for a range of metaUicities, up to 
1000 km/s, although a small number of low metallicity models above 800 km failed to 
converge after several global iterations. A fiducial model with a velocity of 700 km was 
used here, having achieved a converged solution over a range of metaUicities. 

If the shocks are radiative, and if we assume, for example, a density of 1 cm~^, then 
the LycK emissivity will be 0.13 erg cm^^s"^, almost independent of the metallicity of the 
gas. A quarter of this comes from the shock itself, the rest coming from photoionization 
of the assumed neutral gas in front of the shock (the precursor region). If the shock is 
face-on to our sight-line, this would correspond to an observed surface brightness of 1.6 x 
10~^^erg cm^^s^^arcsec^^: nearly two orders of magnitude greater than the observed Lya 
surface brightness (Section 3). Much of this may be self absorbed, but if even 1% escapes. 
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it can explain our observations. The predicted Lya/C IV ratio is also consistent with the 
observations (Francis et al. 1996). 

How does this vary with the assumed density of the gas? If the shocks are radiative, the 
emissivity will scale roughly with the density. Thus even with our lower limit on the density 
p — 10~^cm~^ (Section 5), shocks can produce all the Lyct we observe. 

Would these shocks be radiative? For our lower hmit on the density, 700km s~^ shocks 
become radiative on a timescale of ~ 10^ years, which is comparable to the dynamical 
timescale. Thus it seems hkely that at least some fraction of the shocks will be radiative. 

Thus shocks are probably present, and can easily produce the Lya emission we observe. 
A number of factors could, however, suppress the Lya emission from shocks. 

1. The high optical depth in Lya may have caused us to overestimate the velocity dis- 
persion of the Lya nebula. The emission from shocks is a strong function of their 
velocity: 200km s~^ shocks emit two orders of magnitude less Lya per unit area than 
our fiducial model. 

2. The density of the Lya nebula may be close to our lower limit of p = 10~^cm~^, in 
which case the Lya emission, after the probable self absorption losses, could be an 
order of magnitude below the observed Lya surface brightness. 

3. The Lya nebula could consist of small dense 'bullets' of neutral gas, moving at high 
speeds through a hot, highly ionized, low density medium. The bowshocks in the 
ionized medium would not be radiative, while the reverse shocks within the 'bullets' 
would be slow, and hence have low surface brightnesses (as well as having a small area). 

5.2.2. The Energy Source 

If the shock model is correct, the bulk gas motions explain the Lya luminosity. What 
then could be the energy source for the bulk gas motions? 

If we divide the total kinetic energy of the nebula (assuming a mass of ~ 10^°Mq) by 
the Lya luminosity, we can derive an upper limit on the damping time, which is ~ 10^ years. 
This is less than the dynamical timescale and stellar age of Bl, suggesting that continued 
energy input is required. Note that the Lya we observe may only be a small fraction of 
the Lya emitted, which in turn will only be a small fraction (~ 1%) of the total energy 
dissipation in the shocks. On the other hand, the gas mass may be greater than we are 
assuming. 
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If continued energy input is required, what can the source be? Perhaps the most natural 
energy source is the gravitational potential energy of Bl. If we assume that Bl has a total 
mass of ~ then the gravitational potential energy released in its formation would 

be ~ 10®^ ergs, which is more than sufficient to power the Lya luminosity for the dynamical 
timescale. Such a mass for Bl would imply virial velocities in the Lya nebula of ~ 500km s~^, 
which are comparable to the bulk gas motions observed. This energy could be released by 
the merger of Bla and Bib, or by a continuing mass infall rate of ~ lO^MQyr"^. Haiman, 
Spaans & Quataert (2000) point out that chemically primordial gas at a temperature of 
~ 10^ K has few cooling mechanisms other than Lya emission. Thus if most of the gas 
around Bl is primordial, a Lya luminosity comparable to the gravitational binding energy 
(Section 5.2.2) of the system must be radiated. 

Alternatively, Taniguchi & Shioya (2000) suggest that galactic superwinds driven by a 
starburst could power extended Lya nebulae. This would naturally explain the similarity 
between the dynamical timescale and the inferred age of the stellar population in Bla and 
Bib. The maximum plausible amount of energy which star formation can inject into the 
intergalactic medium is ~ 10^'^crg Mq^ (eg. Bower et al. 2000). As Bla and Bib have a 
combined stellar mass of ~ lO^^Mo, the total energy liberated would be ~ 10 ergs, which 
could in principle drive the observed Lya luminosity for ~ 10^*^ years. 

The Lya emission in high redshift radio galaxies has sometimes been ascribed to mechan- 
ical energy deposition by a radio jet (eg. Bicknell et al. 2000). Could such a jet be present in 
Bl? We do not detect any radio emission from Bl (Francis et al. 1996). Our 3a upper limit 
of 0.23 mJy at 2.4 GHz corresponds to a luminosity limit of lO^^ergs s^^Hz^^ at rest-frame 
1.4 GHz (assuming a typical radio galaxy radio spectrum of the form F^, oc i/~°-^). Typical 
ratios of the total jet mechanical power to the monochromatic radio flux at 1.4 GHz are 
~ 10^^— lO^^Hz-i (Bicknell et al. 1998). Thus we can place an upper limit of ~ lO^'^erg s'^ 
on the energy injection from a radio jet. This is sufficient to power the observed Lya lumi- 
nosity only if > 10% of the mechanical energy is converted into observable Lya emission, 
which seems implausible. Deeper radio observations should test this hypothesis. 



5.3. Other Models 

If the shock emission is suppressed, is there an alternative mechanism capable of pro- 
ducing the Lya emission? 

Stellar ionisation (by hot young stars) is incapable of producing the observed G IV 
emission. Photoionisation by the UV continuum radiation from a QSO can, however, produce 
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both the Lya and C IV, as discussed by Francis, Woodgate & Danks (1997). We do not see 
strong UV continuum emission from any component of Bl, so such an AGN would have to 
be concealed from our sight-line. Chapman et al. (2000) detected strong sub-mm emission 
from a Lya blob at z = 3.09, suggesting that a dusty AGN or starburst was present, but we 
have no evidence for such a source in Bl. 

Extended emission-line nebulosities are associated with many low redshift cooling flow 
clusters, such as NGC 1275. As all high-redshift Lya blobs seem to lie in cluster environ- 
ments, could the Lya luminosity be associated with a cooling flow (eg. Fabian et al. 1986)? 

Steidel et al. (2000) suggest that Lya nebulae are caused by gas cooling from higher 
temperatures, in such a cooling flow. If each cooling hydrogen atom produces one Lya 
photon, however, we require that ~ lO^MQyr"^ of hydrogen be cooling onto Bl. This would 
build the observed baryonic mass of Bl in < 10^ years. This is far smaller than either the 
dynamical time or the age of the stellar population in Bla and Bib, so the simple cooling 
model seems implausible. 

Note that the nebular emission-line luminosities of low redshift cooling-flow clusters are 
also orders of magnitude too great to be explained by this simple recombination mechanism 
(eg. Voit & Donahue 1997). The cause of this discrepancy is still controversial even at low 
redshifts. 

On purely energetic grounds, could the Lya luminosity be powered by the thermal 
energy in the nebula? If wc divide the thermal energy of IO^^Mq of hydrogen at 10^ K by 
the observed Lya luminosity, we get a timescale of ~ 10^ years: shorter than the dynamical 
timescale or the stellar age. lO^^M© may however be an underestimate of the true gaseous 
mass of the Lya nebula, particularly if the neutral gas we observe is immersed in a much 
more massive cloud of hot, highly ionised X-ray emitting gas. 

6. Conclusions 

So what is Bl? On balance, we conclude that Bla and Bib are probably young elliptical 
galaxies, perhaps analogues of E-l-A galaxies today. Both their red colors and compact sizes 
are naturally explained by this model. An extremely dusty starburst in the central regions of 
a disk galaxy can however, with some fine tuning, reproduce the colors and sizes. The bluer 
colors of the centres of both components, and the Lya emission from the core of Bib are, 
however, arguments against the presence of large amounts of dust in these systems. Near-IR 
H-band spectroscopy on an 8m-class telescope should resolve this ambiguity, by detecting 
the spectral signatures of a Balmer/4000A break. Blc appears to be a region of relatively 
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unobscured star formation, perhaps induced by the interaction of Bla and Bib. 

If Bla and Bib are eUiptical galaxies, their surface brightnesses are an order of magni- 
tude greater than those that the Kormendy relation would predict for galaxies of this size, 
and they show strong radial color gradients, with blue cores. The very existence of massive 
elliptical galaxies at this redshift, if confirmed, is interesting (eg. Treu & Stiavelli 1999, and 
references therein). These galaxies would give us a chance to study young elliptical galaxies 
for the first time. Our observations would pose several puzzles: 

• Why is the surface brightness of these galaxies so great? Is this simply due to passive 
evolution? 

• Why do both galaxies show strong, but very similar color gradients? 

• Is the remarkable similarity and proximity of the two galaxies coincidental, or is it 
telling us something about how cluster elliptical galaxies form? 

• If these galaxies are so massive and red at z — 2.38, what would they have looked like 
at z ~ 3 when they were forming? 

The nature of the Lya nebula is less well constrained by our observations. The nebula 
could be excited by photoionization from an AGN, by whatever powers low redshift cooling 
flow nebulae, or by fast shocks powered by either superwinds or gravitational potential 
energy. Spatially resolved spectroscopy of other emission lines will help discriminate between 
these models. 

We wish to thank Joss Bland-Hawthorn and Catherine Drake for their assistance at 
with the TTF observations, and all the staff at CTIO for heroic efforts to get the mosaic 
camera working for us. We are grateful to Mike Dopita, Megan Donahue and Martijn de 
Kool and Ian George for helpful discussions. Hsiao- Wen Ghen was instrumental in getting 
the shtless STIS spectrum extraction software working. David Hogg's cosmology crib sheet 
(astro-ph/9905116) was very helpful to us. 
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